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We propose that the observed matter-antimatter asymmetry can be naturally produced as a 
by-product of axion-driven slow-roll inflation by coupling the axion to by-product neutrinos. We 
assume that grand unihed theory scale right-handed neutrinos are responsible for the masses of the 
standard model neutrinos and that the Higgs is a light field during inflation and develops a Hubble- 
scale root-mean-square value. In this setup, the rolling axion generates a helicity asymmetry in 
by-product neutrinos. Following inflation, this helicity asymmetry becomes equal to a net lepton 
number as the Higgs condensate decays and is partially reprocessed by the SU{2)l sphaleron into 
a net baryon number. 


Axions are attractive candidates for the inflaton be¬ 
cause an approximate shift symmetry protects their po¬ 
tential from large radiative corrections [i|. However, the 
simplest models are difficult to realize within UV com¬ 
plete theories because a Planck-scale axion decay con¬ 
stant is required in order to match the amplitude and 
scale dependence of the observed curvature fluctuations. 
Recently, interest in axionic models has been revived with 
the realization that monodromy effects Q can generate 
a suitable potential for large-field models of inflation. 
These large-field models may be required if any of the 
B-mode signal found in the BICEP2 results Q is due to 
primordial gravitational waves. 

Inflation is efficient at accounting for the adiabatic, 
nearly scale-invariant spectrum of curvature fluctuations 
in the early Universe with an amplitude measured to be 
^ 2.2 X 10“® [3]. In addition, there is an observed 
abundance of matter over antimatter, quantified in terms 
of the baryon-to-photon ratio. Observations of the mi¬ 
crowave background constrain the baryon asymmetry pa¬ 
rameter to be @ 

^ = Ufc - Ufc = ^ 0 15 X iQ-io^ 

and constraints on nucleosynthesis require ry ^ 5.7 —6.7 x 
10“^° in order to get the observed light elemental abun¬ 
dances correct (for a review and references see @). This 
number is remarkably similar to the observed amplitude 
of the dimensionless power spectrum of curvature fluctu¬ 
ations, which leads one to look for a common origin. 

In this Letter we explore the possibility that a net lep¬ 
ton number, sufficient to explain the matter-antimatter 
asymmetry, may be generated via the production of left- 
handed standard model neutrinos during inflation. This 
scenario is naturally accommodated in axion-inflation 
scenarios, involving a dimension-five derivative coupling 
of the axion to by-product neutrinos. In our scenario, 
we assume neutrino masses are generated via the see¬ 
saw mechanism, whereby the standard model neutrinos 
couple to GUT-scale right-handed Majorana neutrinos. 
During inflation, this seesaw may also be active, since 
the electroweak symmetry can be broken by quantum 


vacuum fluctuations during inflation [7|. These quan¬ 
tum fluctuations can generate a Hubble-scale root-mean- 
square (rms) value for the Higgs field which leads to 
the generation of masses for the standard model fields, 
including neutrinos. The axion-inflaton couples deriva¬ 
tively to the neutrino fields, which leads to the asymmet¬ 
ric production of neutrino helicity states, while the see¬ 
saw mechanism ensures that production of right-handed 
neutrinos is highly suppressed compared to the produc¬ 
tion of left-handed neutrinos. Following inflation, the 
Higgs condensate decays, the neutrinos become massless, 
and the resulting helicity asymmetry becomes equivalent 
to a net lepton number. The result is the net production 
of lepton number in the form of left-handed by-product 
neutrinos. The electroweak sphaleron will act on this net 
lepton number L conserving B — L, where B is baryon 
number, while driving B -\- L to zero and thus generating 
a net baryon asymmetry. 

The use of rolling scalars coupled to fermionic currents 
in models for baryogenesis has a long history going back 
to original work by Cohen and Kaplan [^, Q . Dolgov and 
Freese [l^ proposed a rolling axion coupled to the B — L 
current to generate a baryon asymmetry in the presence 
of baryon number violating processes during reheating. 
In contrast to this, our model does not require a complex 
axion that carries baryon number. As we will explain, in 
our model the baryon asymmetry is produced via lepto¬ 
genesis which originates from a helicity asymmetry pro¬ 
duced during inflation. Related ideas have appeared re¬ 
cently such as Higgs relaxation leptogenesis [ill ll2| and 
axion-oscillation leptogenesis [l^ . 

We work in natural units c = h = ks = ^ and denote 
the reduced Planck mass by (SttG). We use 

the 2-component spinor conventions reviewed in 0 , and 
work with a Friedmann-Robertson-Walker metric with 
“mostly plus” convention. 


AXION INFLATION WITH MAJORANA 
FERMIONS 

In addition to the usual action for the standard model 
of particle physics, we consider a model of pseudoscalar 
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inflation together with a set of Majorana fermions 


s = I d^x{v^[^R - - V[4>) 


+ ( 2 ) 


Here (j) is the real pseudoscalar (axion) inflaton with a 
potential V{(j)) which softly breaks the axion shift sym¬ 
metry and drives a period of slow-roll inflation. The Ma¬ 
jorana fermion fields are Vi, which we have rescaled by 
their conformal weight in order to write the deriva¬ 
tive as a partial (rather than covariant) derivative, e^a 
are vierbiens 1^ , while / is a mass scale associated with 
the axion and C is a dimensionless coupling. 

We have written a generic Majorana theory here, but 
we have in mind the neutrino sector of the standard 
model augmented with heavy right-handed neutrinos 
(with Majorana mass terms) to give mass to the by¬ 
product neutrinos via the seesaw mechanism. For sim¬ 
plicity, we will later consider only a single generation of 
neutrinos, and write both the right and left-handed neu¬ 
trino fields as left-handed spinors; = Vi [i < d). 
The physical right handed neutrinos are Npn = v\ 
(3 < * < 6). The axial current is conserved for massless 
fermions, therefore the axion coupling in Eq. m has no 
effect in this limit, contributing only a boundary term to 
the action. 

In local thermal equilibrium, the rolling axion acts as 
a chemical potential for helicity. While the Universe is 
not in thermal equilibrium during inflation, the effect of 
the coupling of the axion to the neutrinos during these 
epochs has a similar effect, and biases the gravitational 
production of one helicity over the other |lfil |. For Dirac 
fermions, where the masses are degenerate rrii = m, con¬ 
servation of charge associated with the additional C/(l) 
symmetry means that this particle production results in 
a helicity asymmetry but not a matter-antimatter asym¬ 
metry. In contrast, for Majorana fermions, with lepton- 
number violating mass terms, this coupling leads to a 
helicity asymmetry which is equivalent to lepton num¬ 
ber and a matter-antimatter asymmetry if these fields 
subsequently become nearly massless. 


these parameters are constant at tree level, they are mod¬ 
ified by both loop and finite temperature corrections. 
The experimentally preferred top quark and Higgs bo¬ 
son masses give loop corrections that result in a neg¬ 
ative running of the coup ling A at a sufficiently large 
vacuum-expectation-value [l7|. For the central values of 
the standard model parameters the electroweak vacuum 
is metastable, however, for values within the 2-sigma re¬ 
gions, a stable potential can be achieved. 

During inflation, in the absence of new physics that 
significantly changes the running of A between the elec¬ 
troweak and inflationary scales, this negative running 
means that the Higgs field is generically light. Conse¬ 
quently, quantum fluctuations of the Higgs field that are 
produced during the inflationary epoch are sufficient to 
generate an rms value 0E1- It is then quite natural 
to assume that the Higgs field generically has a large 
rms value during the inflationary epoch, whose size is 
expected to be of the order 


where A* is the Higgs self-coupling evaluated at the in¬ 
flationary energy scale. The Higgs condensate breaks 
electroweak symmetry during inflation, and the stan¬ 
dard model fields acquire masses set by the Hubble scale 
through the Higgs mechanism. 

For simplicity, we will consider only a single left- 
handed, and a single right-handed neutrino and take the 
form of the neutrino mass matrix to be 


rriij 


0 mo \ yh 

mo M )' “ 72 ’ 


( 5 ) 


where y is the Yukawa coupling and M ~ 10^® GeV is 
the right-handed neutrino mass. 

To proceed, we work in a basis of mass eigenstates, we 
diagonalize the fermion sector by rotating the fermion 
fields. After diagonalization, for a mass hierarchy M ^ 
niD, the seesaw mechanism results in a mass matrix with 
masses rrii ^ M,m\,/M. We assume that the right- 
handed neutrinos are heavy compared to all scales of 
interest, and can be safely neglected. In what follows 
we consider only the left-handed neutrinos, whose mass 
we take to be a free parameter of order the Hubble scale. 


THE STANDARD MODEL HIGGS FIELD 
DURING INFLATION 

The standard model Higgs, $, has a tree-level potential 
of the form 

U($) =A($1'd>)2. (3) 

At the electroweak scale, the parameters /r and A yield 
a stable minimum at a VEV of uew = 246 GeV. While 


LEFT-HANDED NEUTRINO PRODUCTION 
DURING INFLATION 

As noted above, for simplicity we focus on a single 
generation of neutrinos. Varying the action with respect 
to yields the equation of motion for v 

vp = ( 6 ) 
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where m = m'jj/M. We expand each field into a Fourier 
basis, 


where are the Whittaker W functions, r = 

6 and 6' are arbitrary phases, and 


A 




(7) 


where we have introduced creation and annihilation op¬ 
erators, and which satisfy the anticommuta- 

tion relations } = (27r)^d^(k — k')dAA', with all 

other anticommutators vanishing, as usual. We quan¬ 
tize the fields by imposing the anticommutation relations 
{i^q(x, t), 7 r^(y,<)} = iS^S^(x — y). The canonical mo¬ 
menta of the fermions are found in the usual way 


C (j) I 


^ = + (13) 

In order to determine the particle number at any given 
time, we perform a Bogoliubov transformation, and com¬ 
pare the exact wave functions to an instantaneous WKB 
solution that diagonalizes the Hamiltonian. This instan¬ 
taneous quasiparticle number is [illiS 


){k) = 




wx{k\ -I-wa) 


(14) 


TT 


/3 

u 


dC 

dvp 


=iv 


t . ^Od/3 

OL^ ^ 


( 8 ) 


where an overdot here and throughout denotes a deriva¬ 
tive with respect to cosmic time, t. We work in a basis 
of helicity eigenspinors, which satisfy 

a-kx = X^x, A = ±l, e-A(-fc) =iA 6 (fc), (9) 


where the effective frequency and wavenumber are 

ujlit) =kx(tf + , kx{t) = . (15) 

Using the analytic expressions for the wave functions 
from Eq. m, we derive analytic expressions for the 
quasiparticle number during inflation. Assuming that 
m 7 ^ 0 and taking the limit k/aH —> 0, we find 


where lX is a phase that satisfies LX*iX = \ and ■ = 

k k k ' —k 

—lX, Writing the spinors in this helicity basis as 

xkKt) =xkmx{^), y"^“(k,t) = Ytmxi^), (10) 


canonical quantization requires that the fermion wave 
functions satisfy 

^ [x^mkit)+Ytm.kt)] = 1 . ( 11 ) 

A 


Substituting Eqs. dm and © into Eq. ®, the equations 
of motion for the fermion wave functions are 




X^{t) =mYt{t), 


Ytit) =mX^{t). 


( 12 ) 


Note that the effect of the axion coupling is a helicity- 
dependent shift in the effective wavenumber of the modes 
k/aX k/aX + C4>/f- 

Assuming approximate de Sitter space during infla¬ 
tion, H Ri const., and taking (j)/H r; const., the canoni¬ 
cally normalized solutions X^(t), Yu* it ) of Eq. (fl^ that 
match onto the Bunch-Davies vacuum llfill are 


Xkkr) = - — 


im\ 




H ) 


pikr), 


YkikT) = - 


im 

ll 


giS'g A 51? 

VSfcr 




~:{k) =e 




sinh 




(16) 


While Eq. (ITSl) is derived in the limit k/aH —>■ 0, for 
m <C dH it is a good approximation for the particle 
number for modes that satisfy k/aH < d. Note that in 
the absence of the coupling to the axion (d = 0), produc¬ 
tion of both helicity states is symmetric, as expected, and 
highly suppressed for fermions with masses larger than 
the Hubble rate. Eor small masses, the occupation num¬ 
ber approaches its maximum value of 1/2 as m/H —> 0. 
However, for m = 0, the theory is conformally equivalent 
to a Minkowski-space theory, and no particle production 
occurs [ilj. When the coupling to the axion is switched 
on, the particle production is asymmetric between the 
helicity states. For d > 0 (d < 0), particle production of 
the A = -I- (A = —) helicity state is enhanced while parti¬ 
cle production of the A = — (A = -|-) mode is suppressed. 
Larger couplings allow the production of fermions with 
increasingly large mass. 

This particle production can be understood by study¬ 
ing the evolution of Eq. (TT^ in the WKB approximation. 
For m/H <C d, particle production occurs near the points 
of nonadiabatic evolution of k. As discussed in [l^, these 
occur whenever k = 0, where kX/a = —{C/f)4>. After 
this event, k/aH < d, and the resulting (quasi-)particle 
number is approximately constant. The maximum co¬ 
moving wavenumber of excited fermion states can be 
calculated as the maximum comoving wavenumber that 
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makes k = 0. By simple inspection of the terms in fc, we 
can see the scaling fcmax oc {C/f)^. The proportionality 
factor depends on the axion potential, as explained in 

S- 

At the end of inflation the resulting helicity asymmetry 
is (assuming d > 0) 


we used the Friedmann equation to relate the energy den¬ 
sity at the end of inflation to the energy density at re¬ 
heating 




«e \ 

Or/ 


3(l+uj) 




(19) 


n!) = 


E' 

A=±l 


3A 


n^{k)k^dk 


{<) (Ccj^o 


2Tr^a^ \ / 


-QoUf 


(17) 


where is the Hubble rate at reheating, and pe = 
3Mpji7g is the energy density at the end of inflation. Us¬ 
ing Eq. (US, we arrive at an estimate for the asymmetry 
parameter at reheating 


where we have used fcmax oc {C(j)o/f)aeHe, and Oe and Hg 
are the scale factor and Hubble rate, respectively, at the 
end of inflation. The quantit y (n (l) is the phase-space- 
averaged occupation number [l^ of the helicity A. We 
have taken (j)/H = (j)o, where (po ~ Mpi is the field value 
at which inflation ends and oscillations begin. We have 
also summed over the three generations of neutrinos. 

After inflation, once the Hubble rate drops below the 
mass of the Higgs field, the Higgs condensate will de¬ 
cay restoring electroweak symmetry Q and this helicity 
asymmetry will become equal to lepton number. 

Following inflation the axion oscillates which results in 
the production of both helicity states of the left-handed 
neutrinos. However, the helicity states that are produced 
during inflation (and during the first axion zero crossing) 
are produced out to a larger wavenumber due to the fact 
that the other helicity is not produced until the axion ve¬ 
locity changes sign. Provided that the Higgs condensate 
does not decay immediately following inflation, the subse¬ 
quent production events are less efficient due to Hubble 
damping of the axion velocity. The end result is that, 
even after the axion oscillations are taken into account, 
a helicity asymmetry of the order of Eq. (II3 is generic 
for a wide range of parameter space. A detailed analy¬ 
sis of fermion production during and after axion inflation 
appears in |l6| . 


BARYON-TO-PHOTON RATIO 


VR 


(n^) ( Cpo 


Hg 


Tr f ^ 


\ f J Mpi Mpi \ Qg 


Zw 


( 20 ) 


Following reheating, by-product sphaleron processes re¬ 
distribute this asymmetry between lepton and baryon 
numbers [^ . Standard model entropy generation and 
the redistribution of the lepton number into the baryon 
number lowers pp by 1 or 2 orders of magnitude, which 
implies we require ? 7 r ^ 10“^ — 10“® in order to explain 
the present day baryon asymmetry. 

In order to obtain the correct asymmetry, assuming 
that the average occupation number (n+) is of order 
unity, and that the Universe is matter dominated during 
reheating w = 0, we see that we require both high-scale 
inflation and a relatively high reheating temperature in 
order to overcome the suppression. However, the 

tensor-to-scalar ratio constrains the inflationary energy 
scale to be < 10“® — 10“® Mp\. On the other hand, in 
order to prevent the washout of the lepton number stored 
in the left-handed neutrinos, we need to ensure that re¬ 
heating occurs below the scale at which lepton-number 
violating processes involving the exchange of heavy right- 
handed neutrinos are in equilibrium. In order to prevent 
excessive washout, the reheat temperature must be below 
Tr < 3 X GeV 

Taking Tr ~ 10^^ GeV, GUT scale inflation Hg ^ 
10“® Mpi, we see that for / ^ 10“^ —10“^ Mpi we require 
order unity couplings C ^ 1 between the axion and by¬ 
product neutrinos to generate a baryon asymmetry of the 
appropriate size today. 


During reheating, the energy density p = in 

the inflaton is converted into radiation, over an epoch 
where the scale factor expands from ag at the end of in¬ 
flation to Or at the end of reheating. Taking an arbitrary 
equation of state (w) of the Universe between the end of 
inflation and reheating, we write the comoving entropy 
at reheating as 


^rS — -^9*(^R"rR 


Tr 



(18) 


where Tr ^ i/F/Afpi is the reheat temperature, F/ the 
inflaton decay rate, and g* is the effective number of rela¬ 
tivistic degrees of freedom. In obtaining this expression, 


DISCUSSION AND CONCLUSIONS 

In this Letter we have studied the generation of lep¬ 
ton number via the biased production of left-handed neu¬ 
trino helicity states during inflation. In our scenario, we 
have assumed that quantum effects break electroweak 
symmetry during inflation and the by-product neutri¬ 
nos are given masses by GUT-scale right-handed Majo- 
rana neutrinos. Following inflation, electroweak symme¬ 
try is restored as the Higgs condensate decays and the 
left-handed neutrino helicity becomes equivalent to a net 
lepton number. Sphaleron processes redistribute this net 
lepton number, T, into baryon number B, via processes 
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that violate B + L but conserve B — L resulting in a net 
baryon number. 

The number density of neutrinos produced by this 
mechanism is proportional to the cube of the Hubble rate 
at the end of inflation. This is simply due to the nature of 
inflationary particle production, which populates states 
with momenta near the Hubble scale. High-scale inflation 
at or near the GUT scale produces the smallest allowable 
Hubble length and, therefore the largest number density 
of neutrinos. In order that these neutrinos are not diluted 
too much by the subsequent expansion between the end 
of inflation and the onset of the hot big-bang phase, the 
reheat temperature needs to be high. However, this re¬ 
heating temperature cannot be significantly higher than 
^ 10^^ GeV, as scattering processes mediated by heavy 
right-handed neutrinos can wash out the asymmetry. 

To produce the observed baryon asymmetry with this 
mechanism, we require the Hubble rate at the end of 
inflation He ^ 10“®Mpi, a reheat temperature T/j ^ 
10^^ GeV and a derivative coupling between an axionic 
inflaton and a Majorana neutrino sector with strength 
C/f - 102 - 

A number of details of this scenario still need to be 
explored. First, since the fermions are produced during 
inflation, their backreaction could lead to significant con¬ 
tributions to the curvature spectrum. Second, we have 
made only rough estimates of how much of the helicity 
asymmetry eventually becomes a final baryon asymme¬ 
try. This process is dependent on the details of reheating, 
and requires modeling the decay of the Higgs condensate 
and solving Boltzmann equations. Finally, the possibil¬ 
ity of fermion isocurvature perturbations 23| could lead 
to an independent observable that is correlated to the 
baryon asymmetry. We leave these studies for future 
work. 
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